Scaling limits and aging for asymmetric trap models on 
the complete graph and ^-processes 

S. C. Bezerra *t L. R. G. Fontes * § R. J. Gava H P. Mathieu 

V. Gayrard ^$ 



February 21, 2012 



Abstract 

We obtain scaling limit results for asymmetric trap models and their infinite volume 
counterparts, namely asymmetric K processes. Aging results for the latter processes are 
derived therefrom. 
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1 Introduction 

The long time behavior of trap models and related processes with disordered parameters has 
been the theme of several papers in the recent literature. From the inaugurating work of 
Bouchaud [TJ, where the case of the complete graph was shown to exhibit aging, the same 
as well as other cases were analysed. The model on the complete graph was further studied 
in [2] and [3], with different points of view, and considering distinct time scales. And more 
recently, jl] took up the asymmetric case, which is also the model we study here. 

The trap model in the complete graph is sometimes also called REM-like trap model, due 
to its resemblance to a dynamics for the Random Energy Model (REM |5J). Such a dynamics 
for the REM, on the hypercube rather than the complete graph, was studied in [HJ E], where 
aging results comparable to the ones of Bouchaud were derived. See also (SJ El ED]- Trap and 
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trap-like models associated to correlated energy (mean field) spin glasses have been the object 
of more recent work: a dynamics for the p-spin model was studied in [Til [T2] , and results on 
the GREM-like trap model were obtained in [13]. 

Trap models on Z d have also attracted a lot of interest, in connection with aging as well as 
with localization; see [TJ1 [T51 [TBI [T71 Q2] - results on the asymmetric case were obtained recently 
in [221 l2Ql 121]. Analyses on tori were performed in [221 123] . 

In this paper, we revisit the trap model in the complete graph, described briefly below in 
this introduction, and in full in Section |3j Our goal is twofold: 

1. to propose a representation of the model - in terms of trap depth, rather than location - 
for which scaling limits can be derived in a unified manner in different scaling regimes; 

2. and to introduce the infinite volume processes which result from these scaling limits, in 
particular the asymmetric K process. 

Let us now briefly describe the asymmetric trap model in the complete graph with n vertices. 
This is a continuous time Markov chain on the vertices of that graph, whose mean jump time 
at site x is given by r^ _ct , where a e [0, 1] is an asymmetry parameter, and whose transition 
probability from any site x to any site y is proportional to t£, where {r x } are iid positive 
random variables in the domain of attraction of an a-stable law. The random variable T]r a 
may be interpreted as the depth of the trap at site x. The case a = is that of the symmetric 
model. We call the general case where a G [0, 1] the asymmetric model. Let Y n (t) denote the 
site visited at time t. 

This paper is more immediately related to [3] and [I] , so let us briefly outline our results here 
against the background of the ones of those papers. Whereas in the former reference a scaling 
limit was derived for the symmetric model at times of the order of the deepest trap in the 
landscape, and then aging results were derived for a class of two-time correlation functions of 
the limit model at vanishing times, in here we present similar limit results for the asymmetric 
model. Rather than looking at Y n {t) however, we consider Z n {t) = t y ~^, the depth of the 
currently visited trap. As explained below, this is a convenient representation for taking scaling 
limits, not only at times of the order of the deepest trap in the landscape, which we do here 
using this representation (see Theorem 13. ip . obtaining a limiting process which we denote by 
Z, but at shorter time scales as well. We call Z the asymmetric K process, in allusion to the K 
process introduced in [3]. We further derive a scaling limit result for Z at vanishing times (see 
Theorem I3.2p . obtaining a limiting process Z which is self similar of index 1. The latter fact 
may be interpreted as a fuller aging result for Z, involving the dynamics itself, not only a class 
of correlation functions thereof. Other scaling regimes of Z n may be analysed with the same 
approach, with similar results. 

Scaling limits of asymmetric trap models in the complete graph are also the main theme 
of [1]. In that work scaling limits of the clock process are derived in several scaling regimes 
(essentially all of them: from "order 1", where the volume limit is taken first, and then the 
time limit, to the scale where the model is virtually at equilibrium, including scales in between, 
in particular the ones treated here); occurrence of aging and other dynamical phenomena are 
discussed for each regime. 

One reason to consider a representation like Z n , as we do here, rather than the clock process, 
is that, besides the information on the jump times given by the latter process, Z n provides also 
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location information, absent in that process. For, say, correlation functions which depend only 
on jump times (like the II functions discussed on Subsection 13.31 below; see (I3.59ti3.62j) ). the 
clock process is enough. But other ones require location information, and in those cases the 
clock process is no longer enough on its own. We discuss two such examples in Subsection 13.31 
below. 

Z n and Z, as well as their rescaled versions, and Z also, can be described as functions of two 
related subordinators, the second being obtained as the integral of an independent iid family 
of mean 1 exponentials with respect to the first one. Once we obtain the limit of the first 
subordinator in a given scaling regime, a continuity property of the above mentioned function 
implies a limit result for the original process. Section [5] below is devoted to establishing that 
continuity property (see Lemma I2.ip in a somewhat abstract setting, which may turn out to 
be the setting of similar processes of interest. 

In Section [3] we describe our trap models and K processes in more detail and then, applying 
the auxiliary result of Section (2J we derive scaling limit results for them, as anticipated above, 
the one for the trap model in Subsection 13.14 and the one for the K process in Subsection 13.21 
In the closing Subsection 13.31 we discuss the derivation of aging results for three particular 
two-time correlation functions of Z as a corollary to Theorem 13.21 

2 A continuity lemma about a class of trajectories in D 

Let D be the space of cadlag real trajectories on R + = [0, oo) equipped with the J\ Skorohod 
metric (see e.g. |24| Chapter 3, Section 5). Let N* = {1,2, . . .} denote the positive integers. 

Let S, S £ , e > 0, be nonnegative nondecreasing jump functions in D, i.e., suppose that 
there exist (countable) subsets A £ = {xf, i G N*} and A = {xi, i G N*} of 1R + and positive 
number sequences {j x s, i G N*} and {7^, i G N*} such that 

S £ r = Yl 7l f <oo, S r = 7^<oo, r>0. (2.1) 

i: as? e[Q,r] i: 2*e[0,r] 

Consider {T», % G N*}, a family of i.i.d. exponential random variables of mean 1 and let 

i: x\ e[0,r] i: Xi&[0,r] 

'7L , if t G [He _,T e xS ) for some i„, 
Z e t = { l ° l o l » (2.3) 

0, if t£ [T £ x e_,T e x e) for any i, 



and 



z = j 7x l0 , if t G [T Xiq _, T XiQ ) for some i , ^ 
0, if t £ [r ai _, T Xi ) for any i. 



Lemma 2.1 Let S £ , S, Z £ , Z be as above. If S £ —> S, then Z £ — >■ Z in distribution. 
Remark 1 The probability measure in the statement of the above lemma is the one underlying 
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Remark 2 From jQCT, we see toot Z = 5(5, {T u i e N*}) and Z £ = E(S £ , {T h i e N*}), 
where E is the composition underlying the above definitions. Lemma \2.1\ then establishes a 
continuity property of S in its first argument. 

Proof of Lemma 12.11 

We will assume that there exists R' e M + such that \A D [0, R'] \ = oo. Other cases may be 
argued similarly, with simpler arguments. 

Let r _1 be the (right continuous) inverse of T. Let us fix T > 0. Then one readily checks 
that, given 5 > 0, there exists R ^ A, R > R' , such that 

P(r" 1 (T) > R) < 5. (2.5) 

Given rj > 0, we may choose 5' > be such that 

S R+S > -S R <r). (2.6) 

Let us now enumerate A PI [0, R] = {xi,x 2 , ■ ■ ■} such that j Xl > ^ X2 > From the 

hypothesis, there exists m = m(e), with m — > oo as e — > 0, and an enumeration of A 6 D [0, i?] = 
{xf , x§, . . .} such that as £ — > 

sup |x| - Xi\ ) V I m sup |t£ ? - 7J ) -)> 0. (2.7) 

l<j<m / V l<i<m 1 / 

It follows from this and the hypothesis that, given rj > 0, for all small enough e and 1 < k < m 

k k 

i>fc i=l i=l j>fc i>fc 

We now recall that in the J\ topology, functions are close if they are uniformly close in- 
side arbitrary bounded intervals, after allowing small time distortions (for details see e.g. [24] 
Chapter 3, Section 5). 

Now, given k > 1 arbitrary but fixed, independent of e, let {x\, . . . , x^} be an enumeration 
of {xi, . . . ,Xk} such that {xi < ... < Xk}. This leads to an enumeration {xf,...,x|} of 
{x{, . . . , x|} such that for 1 < i < k 

&i ->■ ^ and 7|e ->■ 7^ (2.9) 

(see paragraph of (12. 7p above). At this point we relabel {Tj} so that 7i, . . . , T& are attached to 
xi < . . . < Xfc and commonly to xf, . . . , x|, respectively, which does not change distributions. 
Let and Z^' 6 ^ be the respective versions of Z and Z e with the relabeled {T^}. 

Let us now take a family of temporal distortions (A e ) = (A|) as follows. For 1 < i < k, 
we consider the time intervals Ii = [t~,ti], where ti = T s . and t~ = T s ._, and where 
£f = TL and £?~ = rL_, and then define 



( %r S, if < S < tl 



\%s) 



% i 

p^ts-td+q, ifti<s< tr +v 
( s - *fc+i) + tl+n if s > t fc+1 , 
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(2.10) 



where t k+1 := T R , t e k+1 :— T e R . 

At this point, we have two tasks: the first one is to control the slopes of the functions X £ 
and the second one is to control the sup norm of the difference Z^ e fl — z[ k \ 

We start by the second task. Let M. = U 4 fc =1 /j. If t G M, then 

|4 fc) -^|<mg| 7 , 8 -7l f l, (2.H) 

which goes to zero as e goes to zero by (12.91) . 

If t G [0,tjj~ +1 ] \ Ai, then we have that Z t < j Xk+1 and Z^fl < max^i 7^. Hence, 

I Z[ k) - I < 7 , fc+1 V max 7^ < 7 * fc+1 V E 7^ < E ^ + 2r ?' ^ 

i>k i>k 

for all small enough e, by (12.81) . 

Now, we solve the first problem by considering two cases: 

1) If s G [t~, ti] for some 1 < i < k, then the slope of X £ is given by 



as e -> 0, by (12T9]) . 

2) If s G [tj, iZ_ x ] for some < i < k, where to := 0, then it suffices to prove that 

ti ->• U (2.14) 

tr^t- (2.15) 

as e — > in probability. 

In all cases, the absolute value of the difference of right and left hand sides is bounded above 

by 

m 

E i^i - ^ + E i^i - ^i Ti - ( 2 - 16 ) 

i=l i>m 

The first term vanishes almost surely as e — > by (12.71) and the Law of Large Numbers, 
and, given r] > 0, the expected value of the second term is bounded above by 

E <t + E 7* ^ 2 E 7* + ^ ( 2 - 17 ) 

j>m i>m i>m 

for all small enough e, where use is made of (12.81) in the latter inequality, and (12.14112. 15p follow 
since m — > 00 as e — > and rj is arbitrary. 

To conclude, given < £ < 1,6 > 0, choose T > — log(£/2), and then R satisfying (12.51) . 
and then 5' satisfying (12. 6ft with r\ = £/4, and then such that J^iyk^f^ < C/2- Choosing now 
X s as in (12.101) . we conclude that 

\imsupF{d{Z ik ' £ \Z {k) ) >()<8, (2.18) 

where d is the J\ Skorohod distance on D (see [21] Chapter 3, Section 5). Since 
and Z^ = Z in distribution for all fixed k and £ small enough, the result follows. □ 
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Let us now explain how Lemma [27TI will be used in the sequel. Our aim is to apply it to a case 
where S £ and S are random objects, in fact subordinators, with parameters that are themselves 
random, which we call environment. Both S £ and S, as well as their respective environments, 
will be independent of {Tj}, and the convergence S £ — > S will hold only in distribution: either 
1) the joint distribution of the environment and the subordinators, or 2) the distributions 
of subordinators given the environment, for almost every realization of the environment. In 
both cases, we may use the Skorohod representation theorem (see e.g. [25] Theorem 3.2.2). In 
case 1) we will first explicitly choose a convenient version of the environment, for which the 
distribution of the subordinator, given the environment, converges for almost every realization 
of the environment; with the modified environment, we are effectively in case 2. We can then, 
by Skorohod representation, in both cases, for each choice of the environment, choose versions 
of the subordinators that converge almost surely, and then we are in the setting of Lemma 12.11 
It is clear that the conclusion of the lemma holds for the original subordinator, where the 
distribution referred to in the lemma is the joint distribution of {Tj} and the subordinators 
given the original environment in case 2, and the modified environment in case 1, for almost 
every realization of that environment in each case. In case 1, the result of the lemma will then 
hold for the overall joint distribution of {Tj}, the subordinators given the environment, and the 
environment. 

Establishing the convergence in distribution of the subordinators is done by verifying the 
convergence of the respective Laplace exponents. 



3 Application to trap models on the complete graph and 
K processes 

We will apply the lemma above to show scaling limit results for trap models in the complete 
graph and for K-processes. We introduce these two processes next. 
We first consider the trap model on the complete graph 

K n = {{l,...,n}, {(x,y), x,y = l,...,n}} (3.1) 

with n vertices (differently from the usual definition, here we include self loops, for convenience - 
this should not matter in the convergence results below): Y n = {Y n (t)) t >o, which is a continuous 
time Markov chain with jump rate at site x given by 



T 



-(1-tt) 



x 



(3.2) 



and transition probability from site x to site y given by 



— , (3-3) 

where a G [0,1] is a parameter, and r := {r x ,x G K n } is an independent family of positive 
random variables with common distribution in the domain of attraction of a stable law of degree 
< a < 1, that is, 

p ( r i > = ^r, * > o. (3.4) 
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where L is a slowly varying function at infinity. 

We call Y n an asymmetric or weighted trap model on the complete graph with asymmetry 
parameter a, mean jump time parameters {T^~ a , x — 1, . . . , n} and weights {r", x — 1, . . . , n}. 
The latter set of parameters may indeed be seen as unnormalized weights of the transition 
probabilities of Y n . Notice that the a = (symmetric) case corresponds to uniform weights. 

We will consider the following construction of Y n . Let 

M = {N^ := (N^ x) ) r >o, x G N*} (3.5) 

be a family of independent Poisson counting processes such that the rate of is r". Let erj 
the j-th event time of N^ x \ j > 1. Let also 

T = {T ; t!> x \ x G N*} (3.6) 

be independent mean 1 exponential random variables, independent of Af and r, and define for 
r > 

n n N$. x) 

x=l x=l i=l 

Then 

Y n {t) = x, if T n (of } -) < t < T n (af ] ) for some xj > 1. (3.8) 

is a construction of Y n as above described, with initial state distributed on {1, . . . ,n} in such 
a way that site x has probability weight proportional to r", x G {1, . . . , n}. 

Remark 3 Regarding the latter point, notice that the initial state ofY n is the one whose Pois- 
son mark is the earliest, so it corresponds to the minimum of n independent exponential random 
interarrival times with rates t£, x G {1, . . . ,n}, and it is well known that the probability that 
the minimum of n independent exponential random variables is a given such random variable 
is proportional to its rate. 

We are interested below in 

z n(t) = 4:; y (3.9) 

This is the representation for the process aluded to at the introduction above. It has been 
considered in [2], where the symmetric (a = 0) case was studied, and a (single time) scaling 
limit result was derived for it, first taking the volume, and the time, to infinity (see Proposition 
2.10 in that reference) - this is an aging regime not considered in this paper, but rather in jl]. 

Remark 4 Z n and Y n may be seen as processes in random environment, where t is the set of 
random parameters acting as environment. Indeed, given r, both are Markovian (this should 
be clear for Y n , but a moment's thought reveals that it is true for Z n as well, even when there 
are same values for t,i 's with distinct i's). Notice also that t is an environment for S n as well, 
which for each n > 1 is a subordinator for every fixed such environment (recall the discussion 
at the end of Section \2.1\) . This aspect, which is characteristic of the complete graph, makes our 
approach particularly suitable, since by an application of (the continuity) Lemma \2.1\ we are 
left with establishing convergence of subordinator s (in the Skorohod topology), which reduces to 
showing convergence of Laplace exponents (in the topology of real numbers), which is relatively 
simple, as we will see below. 
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Remark 5 Given S n , Z n may be identified in distribution to E(S n ,{Ti, i G N*}), with S 
introduced in Remark^ 

We now turn to K-processes, which is a Markov process in continuous time on N* = 
{1,2,..., 00} constructed in a similar way as Y n was above, as follows. Let 7 = {j x , x G [0, 00)} 
be the increments of an a-stable subordinator in [0, 00) given by a Poisson process V in 
(0, 00) x (0, 00) with intensity measure 

ai~ 1_a dx dy. (3. 10) 

Let 

M = {N {x) := (N!; x) ) r >o, x G N*} (3.11) 

be a family of independent Poisson counting processes such that the rate of is 7". Let 

the j-th event time of N( x \ j > 1. Let T be as above (see (13.61) ). and independent of AT. 
Define for r > 

S{r)= £ llr a Ni x \ r(r)= £ 7^ a E T ^' ( 3 ' 12 ) 
xe[o,i] ie[o,i] «=i 



and then make 



1 if r(<3"f -) < i < r(oj* } ) for some x,j > 1, 

00, otherwise. 



Remark 6 It can be verified that when a > a, then Y is a jump process, and so there is almost 
surely no t for which Y(t) = 00 (since in this case Uj iX \T(aj — ), r(aj x ^)) = [0, 00)). And in the 
case where a < a, there almost surely exist t's for which Y(t) = 00. (One way to check these 
claims is by verifying that when a > a, {a^; j > l,x G [0,1]} is a discrete subset of [0, 00) 
almost surely, and when a < a, it is almost surely dense in [0,oo) ; and these in turn follow 
from the fact that J2xe[o 1] ^ s a ^ mos ^ surely finite in the former case, and infinite in the latter 
one.) 

Let 

Z t = ri; a , (3.14) 

where 700 should be interpreted as 0. 

Remark 7 Z and Y may be seen as processes in random environment, where 7 (more specifi- 
cally, 7|[o.i] = {7z? x G [0,1]},) is the environment. Indeed, given 7, both are Markovian. 7 1 [0,11 
is also an environment for S, which is a subordinator for every fixed such environment (recall 
the discussion at the end of Section \2.1\) . 

Remark 8 Given S , Z may be identified in distribution to S(S I , {T,, i G N*}) ; with H introduced 
in Remark^ 
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Remark 9 In and other references the representations used for the trap model and K process 
are the ones given here by Y n (t) and Y(t), t > 0, respectively (see k3. 8\) and h3.13\) above). The 
alternative representation Z n (t) and Z(t), t > 0, we adopt here (see A3.9\) and ft3.14\ ) above) has 
the advantage of leading to a unifying approach for taking the scaling limits of those processes, 
as explained in the introduction and will be done in detail in Subsections \3.1\ and \3.^\ below. 

In the next subsection, we will consider a particular scaling regime for Z n and establish 
a scaling limit result under which Z n converges to the K process. Then, in the following 
subsection we will derive a scaling limit result satisfied by Z. All proofs will rely on Lemma I2TT1 
above to get the results from the convergence of the appropriate S e in each case (see statement 
of that lemma and its preliminaries above). In order to obtain the latter convergence, since 
we have subordinators in all cases, it will suffice to establish convergence of the associated 
Laplace exponents. The last subsection is devoted to a discussion on aging results (for two- 
time correlation functions) satisfied by Z as a consequence of Theorem 13.21 and other results. 

3.1 Scaling limit for Z n at large times 

For r > 0, let 

U(r) = J2 7- (3-15) 

x£[0,r] 

Given a sequence (c n ) n >i, set 

Z[ n) = <i-Z n (t/£-), t > 0. (3.16) 
Theorem 3.1 There exists a deterministic sequence (c n ) n >i such that 

(Zi%> -> {Z t ) t > . (3.17) 

in distribution as n — > oo. 



Remark 10 The distribution referred to above is the joint one of the dynamical and envi- 
ronmental random variables. The sequence (c n ) will be exhibited explicitly in the proof below 
fsee\3~2l\). 

Proof 

By Lemma I2.1[ and recalling the discussion at the end of Section [21 it is enough to establish 
the limit 

S (n) -)• S, (3.18) 

where 

n 

Sf ) := <*-SMr) = X)^,) 1 -^, (3.19) 

2 = 1 

since, given S<> n \ is identically distributed with Z(S( n \ {T u i e N*}). (See the discussion 
at the end of Section [2] above.) 
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Since S^ n \ n > 1, and S are subordinators (given r), it is enough to show the convergence 

of the corresponding Laplace exponents (see Corollary 3.6 page 374 in [26]). A straightforward 
computation yields 

n 

<p n (\) := ]T(c n r x ) a (l - e-^) 1 -"), A > 0, (3.20) 

x=l 

as the Laplace exponent of S^ n \ In order to establish its limit (a weak one), we will make a 
precise choice of c n and switch to another version of r, so that the limit will be strong. We 
follow [H], Section 3. Let 

c n = (inf{t > : P(n > t) < r^ 1 }) -1 , (3.21) 

ri n) := c" 1 g n {U(x) - U(x - l/n)) , x e (0, 1] n (3.22) 

g n (y) = c n G- 1 (n 1 / a y),y>0, (3.23) 

where G _1 is the inverse of the function G defined by the following condition. 

¥(U(1) > G{x)) = P(n > x), x > (3.24) 

Then 

y3 n (A):= (T^ni-e-^) 1 -) (3.25) 
!ce(o,i]niz 

is identically distributed with </? n (A) (simultaneously for all A). Now let 

X 5 = { x g [0, 1] : j(x) >5} = {xi< ...< x K }, (3.26) 

and 

~(n) _ j An) _ 1 r „ i . . (n) I 



Tl j = <^ X™ = - \nx{\ < . . . < x " j = - \nx K ] . (3.27) 
^ n n J 

where the strict inequalities in (I3.27P hold provided n is large enough (for each fixed 5). 
Now Lemma 3.1 in |14j implies that for every 5 > 

£ (r|")r(l - e"^) 1 -) -> £ 7 «(1 - e-^-) (3.28) 

almost surely as n — > oo. One also readily checks that 

£ (ri«)) a (l - e-^n < A ^. (3.29) 

se(o,i]n±zv4 n) xe(o,i]n^\^ n) 

Since, as argued in paragraphs of (3.25-3.28) in [T4] , we have that the lim^o limsup n _ >00 of the 
sum in the right hand side of ( 13 .29 j) vanishes almost surely, we may conclude that 

(p n {\) -> - e" A7 " a ), A > 0, (3.30) 

xe[o,i] 

almost surely. The right hand side is the Laplace exponent of S, so the proof is complete. 
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3.2 Scaling limit of Z at small times 

In this subsection, we assume < a < a. Let 

zf = e- l Z £t . (3.31) 

Before stating a convergence result for Z^ e \ let us describe the limit process. Let (S t )t>o be 
an d-stable subordinator, where 

a — a . 

a = , (3.32) 

1 — a 

and whose Laplace exponent is given by <p(X) = c\ a , where c is a constant to be determined 
below. 

We may then write S as a partial sum of its increments as follows. 

S r = 7* > (3-33) 

a;6[0,r] 

where {%, x G N*} are the increments of S. 
Let now 

xG[0,r] 

where {T x , x G [0,cxd)} is an iid family of mean 1 exponential random variables, independent 
of S. 

Remark 11 One may readily check that T is also an a-stable subordinator (under the joint 
distribution of S and {T x , x G [0,oo)},). 



Now define 



7 _ ■ %, if t G [f x _, F x ) for some x G [0, 00) ^ 
0, for all other t > 0, if any. 



Remark 12 Z may be seen as a process in random environment, where S is the environment. 
Indeed, given S, Z is Markovian. And the distribution of Z (integrated over the environment) 
makes it a self similar process of index 1, that is, (Z t ) t > Q = (c _1 Z ct ) t > in distribution for every 
constant c > 0. This latter property explains the aging behavior of Z in its small time scaling 
regime, as established below. 

Remark 13 Given S, Z may be identified in distribution to E(S, {T i; i G N*}), with E intro- 
duced in Remark^ 

We are now ready to state this subsection's result. 

Theorem 3.2 If < a < a then for almost every 7 

(Zj: £) ) t > -> (Z t ) t > . (3.36) 

in distribution as e — > 0. 
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Remark 14 The distributions referred to above are, to the left, the conditional distribution of 
(Z^) given 7, and to the right, the (unconditional) distribution of (Z t ). 

Remark 15 An argument for the integrated version of Theorem \3.2\ (where the stated weak 
convergence is with respect to the joint distribution of Z given 7 and 7 J is simpler than the one 
given next. See Remarking below. 

Proof 

Let 

^ = e ~' E 7^~m> r > (3.37) 
xe[o,i] 

where a was introduced in (I3.32p above. Then, given 7 ande > 0, (S ( t £ \ t > 0) is a subordinator, 
and its Laplace exponent equals 

^ (£) (A) = e & J2 7, a (l - e"^"), A > 0. (3.38) 
xe[o,i] 

By Lemma HH1 and recalling the discussion at the end of Section [2J to get the result, it is 
enough to establish the limit 

S {£) ->• S (3.39) 

in distribution as e — > for a.e. 7. Since we are dealing with subordinators, it suffices to show 
that for almost every 7 

{e \\) ->• cX & , A > 0, (3.40) 

as e — > 0, for some positive constant finite c. This is obvious for A = 0, so let us fix A > 0, and 
write 

A-V £ )(A) = R~ a (R^Ti 1 ~ e- {R ^ a ) (3.41) 

x6[0,l] 

with R = (e^A) 1 ^, and then argue in the sequel that the left hand side converges to a constant 
as R — >■ 00 for a.e. 7. 

We start by considering 

R5- 1 

W ■= R ~° E E (^) a (! - ^^"X^^-i),^}- (3-42) 
xe[o,i] »=i 

Since the difference between W and the left hand side of f)3.4ip is bounded above by 

R- {a ~ a) E tfWh (3-43) 
xe[o,i] 

which vanishes as R — > 00 for a.e. 7, it is enough to establish the convergence result for W. We 
estimate it as follows. 

RS- 1 R5- 1 



w-x x < R- a x t ■= R~ a E ~ e ~ (Si)1 a ^ Mi ( 3 - 44 ) 

i=2 i=2 
RS- 1 R6- 1 

W > J R- Q ^Xr : = J R-^(5(z-l)) a (l-e- (5(i - 1))1_a )M i , (3.45) 



i=2 i=2 
12 



where X x — R "Eiejoi]^)^ ~~ e a )I{ 7a . e [o,^]} an d Mi is the number of points of V 

in the region [0, 1] x [-|(z — 1), (recall paragraph of (13.101) above). 

X\ can be bounded above by R~ a XLe[o i] (-^7*) ^{ 7:E e[o -]}> an< ^ this nas ^ ne same distribution 
as -R _a XLe[o R a ] 7^ ^{7o:G[o,5]} for every R > 0, by the scale invariance of 7. We can use standard 
large deviation estimates for the latter expression to conclude that X\ can be ignored in the 
limits as R — )■ 00 and then 8 — > (here we may use the existence of a positive exponential 
moment for XLe[o 1] 7^ I^elo^]} f° r an y <5> a result that follows as an application of Campbell 
Theorem - see [27J). We concentrate on the right hand sides of (13.441 13.451) . 

We start with (I3.44p . By the exponential Markov inequality, we get, for given > 0, 

R8- 1 R&- 1 \ . 



P J R- Q ^X+>i?- a ^EX+ + H <- (3.46) 

\ i=2 i=2 J 

where A = Ee 6 ^ 1 x t and B = e «ES~ l '"tf+« a C. 

Since Mj, i > 2, are independent Poisson random variables, we obtain 

— = e -^«0+Efi 2 ~V* 9 -i-^)EM^ 
B 

where a = (5i) a (l - e" (<5i)1 ~ a ). 

We choose 6 = R~ b with a < b < a < 2b. Then, using the estimate 

& a , R a 



EMi = / ——dx < — — - , (3.48) 



! id-i)X 1+a ~ <^(i-l) 
we find that the sum in the exponent in (I3.47P is bounded above by 



m 1 pa oa-26 RS 1 



i=2 V ' i=l 

Since the sum on the right of (13.491) is bounded by constant times R 2a ~ a V log R, and using the 
above estimates, we find that the exponent in (I3.47P is bounded above by 

- R a - b £ + const R~ c \ (3.50) 

for some constant d > 0. We can then apply Borel-Cantelli and conclude that for a.e. 7, given 

RS- 1 R5' 1 

r-° x t < R~ a EX f + £ ( 3 - 51 ) 

i=2 i=2 

for all large enough R. 

Conversely, we can conclude that given £ > 0, for a.e. 7 and all R large enough 

R5- 1 RS- 1 

R- a x r ^ R ~ a E EX * r - ( 3 - 52 ) 

i=2 i=2 
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(13311 and fl332j) then imply that 

R5- 1 R6- 1 

lim inf FT a V EX~ < lim inf W < lim sup W < lim sup iT Q V EX+ (3.53) 

H->00 ^—^ i?-S>00 fl^Qo JJ^oo ^— ' 

To conclude, it is enough to verify that 

lim inf lim inf R~ a V EXr = lim sup lim sup R~ a V EX+ (3.54) 

i=l i=l 



is a (positive finite) constant c. 

We begin with the following estimate. 



EX+ = (<5z) a (l - e -W x - a ) / -^dx < (5i) a (l - e'^ 1 '")- A/ r_ (3.55) 

Summing up: 

Z. - M ^ ( , ( . _ 1))1+a 

- « £ r-°—^r- (j^t) (3*) 

i=2 v 7 

Now as i? —t- 00, the latter sum converges to a series, which is readily seen to be an approxi- 
mation to an integral. We find that 

lim sup lim sup /T a V EX+ < a / — — - da. (3.57) 

8^0 i?^oo . 2 Jo X a a 

We similarly find the latter expression as a lower bound for 

R6- 

nf /T Q 

5-s-O i?^oo 



lim inf lim inf iT Q V EXr 



i=2 



and ( 13.541) follows, with the right hand side of (13. 57ft as the constant c. □ 

Remark 16 In order to get an integrated result, rather than the quenched one above, it is 
enough to replace {7^, x G [0, 1]} by a suitable version. In view of the right hand side of Jff.^ip , 
we consider {R~ lr j xi x G [0,i? a ]} ; and that expression thus becomes 

R- a J2 r x (l~e-^ a ). (3.58) 

xe[o,R a ] 



l-a . 



Since Ylxelo 1] 7"(1 — e ~ Jx ) ^ s integrable, as can be checked by an application of Campbell 
Theorem, the Law of Large Numbers yields a quenched result for the version, which implies the 
integrated result for the original 7. 
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Remark 17 A few words about the cases where a < a < 1. When a > a, we have that Z 
is a jump process in N* (see Remark \@\) with Z(0) = ^ x with probability proportional to 7", 
x G [0, 1]. It follows then that Z^ — > 00 identically almost surely as e — > 0. 

The case a = a demands more delicate analysis. We have that (p^(X) (see \3. 38\) , when scaled 
with a factor of \ loge\~ 1 (instead of e a = 1 in this case), converges to a number r independent 
of A > as n — )■ 00 in probability, and this is the Laplace exponent of a subordinator which 
equals for an exponentially distributed amount of time of rate r, and then jumps to 00, where 
it stays. One may then argue from this that Z^ — > 00 identically as e — >■ in probability. 

3.3 Aging in the K process 

Theorem 13.21 may be viewed as an aging result for Z, since Z is nontrivial and self similar with 
index 1. Corresponding aging results for two-time correlation functions follow. 

Below we consider three examples of correlation functions related to aging, and derive scaling 
limit/aging results for them as a consequence of Theorem 13.21 (as well as of other results derived 
above). Other correlation functions can be similarly treated. 

Example 1 We start with the time correlation function introduced in pQ, which is the one 
that is usually studied in connection with his model. Let 

n(t, s; 7) = P(no jump of Z on [t, t + s] (7) (3.59) 

(see Remark [T8l below) . 

Let $ G D, let £>($) denote the set of discontinuities of $, that is, X>($) = {t > : $(t) ^ 
$(£—)}, and consider F:Dx(0,oo)x (0, 00) — > {0, 1} such that 

F($; t, s) = l{[t, t + s]C] £>($) = 0}. (3.60) 

Then we have that 

n(£t,£s; 7 ) = E[F(Z (e) ;t,s)| 7 ]. (3.61) 

Let also 

fl(t,s) =P(no jump oflon [t,t + s]). (3.62) 

Since deterministic single times are almost surely continuity points of Z, we have that 
F(-;t,s) is almost surely continuous under the distribution of Z. We thus conclude from 
Theorem 13.21 that if < a < a, then for almost every 7 

lim n(et, es; 7) = E[F(Z; t, s)} = fl{t, s). (3.63) 

The aging phenomenon, namely !!(•,•) being a (nontrivial) function of the ratio of its argu- 
ments, then follows from the self similarity with index 1 (and nontriviality) of Z, but in this 
case there is an explicit expression for IT, obtained as follows. One readily checks that the right 
hand side of fl3T62|) equals P([t, t + s] nK(f) = 0), where K($) is the range of $ G D. Since f is 
an d-stable subordinator (see Remark [TTT) . an application of the Dynkin and Lamperti arcsine 
law theorem for that probability yields 

n( t) 8 ) = S J^™1 f e - & (i - ef- 1 de. (3.64) 

71 Js/(t+s) 
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The limit in (I3.63j) was first obtained in [T] (as the expression in ( I3.64p ) for the case where 
a = 0. The general case < a < 1 was first studied in [I] (see Theorem 3.3 for the case a < a, 
and Theorem 3.4 for the case a > a; the particular limit (I3.63j) is (7.5) in that reference). 

In case a > a, then the discussion in Remark [T7] indicates that the limit in ( I3.63P is 
identically 1, and that aging is thus interrupted. 

For the next examples, we restrict a to [0, a). 

Example 2 Let 

R(t, s; 7) = ¥(Z(t) = Z(t + s)| 7 ). (3.65) 
Then, the difference between R(et 1 es\ r y) = P(Z t = Z^l^y) and Tl(et,es; / -f) is given by 

nz[ £) = Z$ B ; ? Z[% for some r E [0, s] (3.66) 

Let and V denote the point processes in (0, oo) x (0, oo) associated to and S, respec- 
tively, i.e., 

= {{t, sl £) - sty-, t > o, sl £) - > o} , v = [(t, s t - s t _y. t>o,s t - s t - > o} . 

(3.67) 

The convergence in distribution S^ — > S argued in the proof of Theorem 13.21 implies that 

V (£) -> V (3.68) 

as e — > in distribution (in the point process sense; for almost every 7). 

Let also T ie \t) = e^T^t), t > (see paragraph of (13712]) above). We have that 

f (e) -> f (3.69) 

in distribution for almost every 7 (see (I3.34p above). This claim may be argued as follows. 
Since (f| ) is a subordinator, an entirely similar reasoning to the one employed in the proof 
of Theorem 13.21 may also be employed to establish this result. It also follows from a continuity 
property of (f^) as a function of (S^) and T similar to the one established in Lemma [2.11 
and similarly proven. We leave the details for the interested reader. 
For arbitrary 5, T > 0, consider now the event 

4ks = {Zt £) > 8, Zit > 5, f( £ )(T) > t + s}, (3.70) 

and let B$j, be the event that there exist two points in R {(0,2T) x (5/2, 00)} with the 
same second coordinate. Now one readily gets from the above convergence results that 

limP(4^ t J 7 ) = F(Z t > 5, Z t+S > 5, f(T) >t + s), (3.71) 

£—5-0 ' ' ' 

and this can be made arbitrarily close to 1 by choosing 5 and T appropriately. We also have 
that 

limP(Bg| 7 )=P(B*,T), (3-72) 
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where B$,t is the event corresponding to Bg £ j, upon replacing by V. The latter probability 
clearly vanishes. Now, since the intersection of the event in the probability in (13 .66j) and A s e ' T t s 
is contained in Bgj,, we conclude from the above that 

lim R(et, es; 7) = lim ftfet, es; 7) = tl(t, s) (3.73) 

£->0 £-S>0 

for almost every 7. 

Remark 18 The aging correlation functions 

n(t,s;7) = P(no jump ofY on [t,t + s]\^), (3.74) 
R(t,s;j) = F(Y(t)=Y(t + s)\-f) (3.75) 

are more widely considered in the literature than their barred versions A3. 59\) and A3. 65\) above. 
In the present case there is almost surely no difference, since a.s. 7 X 7^ j y provided x 7^ y and 
lx > 0. 

The above examples could be done either by considering the clock processes and T on 
their own, together with (I3.69p . in the case of Example 1, or, in the case of Example 2, we used, 
besides Theorem 13. 2} convergence results for S and V (in the appropriate scale), and in both 
examples the limit is a correlation function of the limiting clock process T. Our last example 
is natural from the aging point of view, requires Theorem 13.21 alone, and the limit is not a 
function of F alone. 

Example 3 Let 

Q(t, s; 7) = P (sup re[M Z(r) < sup re[0it+s] Z(r) \j) . (3.76) 

This function was suggested in [H] as a "measure of the prospects for novelty in the system". 
Z is almost surely continuous in single deterministic times, so we have that 

lim Q(et,es;j) = P (sup re[0)t] Z(r) < sup r , e[0jt+s] Z(r)j =: Q(t,s), (3.77) 

since the function l{sup r€ [ t ] $(r) < sup rg [ t+s ] $(r)} is continuous in $ 6 D for almost every 
$ under the distribution of Z. We note that Q(t, s) is a function of the ratio t/s only, by the 
self similarity of Z, but an explicit expression is not available, as far as we know, as it is for 

h(t,s). 
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